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Summary 

The radar scattering properties of realistic aggregate snowflakes have been calculated using the Rayleigh- 
Gans theory. We find that the effect of the snowflake geometry on the scattering may be described in terms of 
a single universal function, which depends only on the overall shape of the aggregate and not the geometry or 
size of the pristine ice crystals which compose the flake. This function is well approximated by a simple analytic 
expression at small sizes; for larger snowflakes we fit a curve to our numerical data. We then demonstrate how 
this allows a characteristic snowflake radius to be derived from dual-wavelength radar measurements without 
knowledge of the pristine crystal size or habit, while at the same time showing that this detail is crucial to 
using such data to estimate ice water content. We also show that the 'effective radius', characterising the ratio 
of particle volume to projected area, cannot be inferred from dual-wavelength radar data for aggregates. Finally, 
we consider the en'ors involved in approximating snowflakes by 'air-ice spheres', and show that for small enough 
aggregates the predicted dual wavelength ratio typically agrees to within a few percent, provided some care is 
taken in choosing the radius of the sphere and the dielectric constant of the air-ice mixture; at larger sizes the 
radar becomes more sensitive to pailicle shape, and the errors associated with the sphere model are found to 
increase accordingly. 
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1. Introduction 

Aggregation plays an important role in the evolution of ice clouds and the devel- 
opment of precipitation. Understanding the geometry and size distribution of ice ag- 
gregates is crucial to analysing radar returns and interpreting them in terms of particle 
size and ice water content. A recent study (Westbrook et al 2004a,b) presented a new 
theoretical model of ice crystal aggregation which produced 'synthetic' snowflakes with 
realistic geometry and distribution by size; sample aggregates from those simulations are 
shown in figure 1 , alongside images of real aggregates from a cirrus cloud. This model 
includes no a priori assumptions about the relationship between snowflake mass and 
dimension, or about the distribution by size. The results from this study show that many 
of the features of the resulting aggregate snowflakes are universal - that is to say that 
they do not dependend on the shape or initial size distribution of the pristine ice crystals 
which make up the flakes. 

The purpose of this paper is to calculate the radar scattering properties of these 
synthetic aggregates and to consider how the results may be used to interpret dual 
wavelength radar data in terms of average particle size and ice water content. We 
identify a number of features in the scattering behaviour of the snowflakes which are 
universal and, equally importantly, a number of features which are not. These results 
have significant implications for the estimation of the microphysical properties of ice 
clouds from radar data. We also consider the errors involved in modelling the aggregates 
as spherical ice/air mixtures, and in particular what size sphere should be used in place 
of the aggregate it is intended to represent. 

* Corresponding author: National Centre for Atmospheric Research, Boulder, Colorado, USA; e-mail 
prfield@ucaredu 
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Figure 1. Aggregates of bullet-rosette crystal types: top half shows images obtained from an aircraft flight 
through a cimis cloud at temperatures between — 44°C and — 47°C (~ 9km altitude), using a cloud particle 
imager. Lower half shows projected images of some of the 'synthetic' aggregates produced by our computer 

simulations (see text). 



Here we present a brief overview of the Westbrook et al (2004a,b) model, and the 
key results from it. The model is based on the rate of hydrodynamic capture between 
two spheres, well known in the raindrop coalescence literature: 



where D is the maximum particle dimension and v is the fall speed. Unlike raindrops, 
snowflakes are clearly not spherical: in order to accurately sample collisions between 
the complicated aggregate shapes we use Fjj as a rate of 'close approach' between the 
two particles i and j. We pick pairs of particles at random with a probability proportional 
to Vij, so as to choose two that are likely to undergo a collision. We then track the pair 
along one of the possible trajectories that the collision area f + DjY encompasses. 
If a collision occurs, the particles are stuck together rigidly at the point of initial contact; 
if the pair miss one another, they are returned to the system and a new pair is picked. In 
this way we correctly sample the collision rates between the complex aggregates. 

The model is completed by an explicit form for the particle fall speeds, and this is 
provided by Mitchell (1996): 



where mg is the particle weight and r is a characteristic radius (see below). The 
parameter a governs the hydrodynamic regime: here we assume an approximately 
inertial flow and set a = ^. The air enters the expression through its density p and its 
kinematic viscosity Vk- The characteristic radius in equation|2lis the radius of gyration 
(see appendix A); however systematically replacing r with the maximum dimension D 
(as per Mitchell) does not affect the results described here. Our computer simulations 
show that r (X D, and since we are only interested in how likely one collision is relative 
to another (in order to pick pairs of particles), it is unimportant which characteristic 
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Figure 2. Snowflake mass as a function of radius of gyration for synthetic aggregates of 50 and 100/^m columns. 
Mass is in units of mo (the average mass of the pristine columns). The simulation data (circles) quickly approaches 

the asymptotic m <x scaling (dashed line). 

length scale we use to calculate v. From the point of view of the radar scattering results 
the radius of gyration is a particularly natural length scale, and in what follows we will 
in general use r rather than D. For reference, our simulation results show that r ~ 0.31?. 

3. Snowflake geometry and size distribution 

The aggregate snowflakes produced by the model have a fractal geometry (in a 
statistical sense). Consequently the flakes have a rather open structure, with a power law 
relationship between mass and radius: 

m = ar'^f. (3) 

where the exponent dj is the fractal dimension and has value of less than three. This 
kind of power law scaling has been widely reported in the literature, usually with 
an exponent of around two (eg. Heymsfield et al 2002: = 2.04 for aggregates of 
bullet-rosettes, df = 2.08 for aggregates of side -planes; Locatelli and Hobbs 1974: 
df = 1.9 for aggregates of plates, side-planes, bullets and columns; Mitchell 1996: 
df = 2.1 for aggregates of side-planes, columns and bullets). This tallies well with our 
simulation results where we measure df = 2.05 it 0.1, and our theoretical arguments 
(see Westbrook et al 2004b) which lead to the prediction df = \/{l — a) = 2. 

Crucially this asymptotic scaling is determined purely by the hydrodynamic regime 
and not the initial conditions (pristine crystal type or size). These details do have an 
affect on how quickly the asymptotic regime is approached, but as figure|2lillustrates, in 
general only a few collisions are needed before equation|3lapplies. 

Unlike df, the prefactor a is not universal, and contains all the information relating 

to the (average) mass mg and radius tq of the pristine particles: a = 'jmo/rQ^ . The 
dimensionless geometrical factor 7 is determined by the crystal shape. 

The size distribution for our synthetic aggregates was found to approach a universal 
underlying shape, which spreads out as the aggregation continues: 

n(m, t) = m~j(j){m/mav) (4) 
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Figure 3. Size distributions by mass for our synthetic snowflakes. Left panel shows snapshot distributions at 
three different times in the evolution of the size distribution: niav {t) = 20 (black line), 55 (dark grey line) and 150 
(light grey line). Mass is in units of mo as per figure|2| the simulation started with 5000 columns of unit mass, 
and 5000 with twice that weight. In the right hand panel we plot m|„n(m, t) against rajmav to show that all 
three are in fact the same underlying distribution rescaled as a function of mav For our simulations total mass 

is conserved and as a result we take 5 = 2. 



where we define the weight average mass mav{t) = Yl I ™^ ^s a positive 
constant. The function ^ describes the underlying distribution shape and its form is in 
principle sensitive only to the value of a in the velocity law As the aggregation 
progresses, the average mass increases, stretching the distribution whilst reducing 
the overall concentration through the factor ■ Equation |3 reflects the idea that 
there is a single underlying distribution which is rescaled depending on how far the 
aggregation has evolved (as characterised by the average aggregate mass viiav)- 

To test this 'dynamical scaling' we measured n(m, €) at three different points in the 

evolution of our simulations, and plot 'mivn{m, t) as a function of m/mav in each case. 
The results are shown in figure|51and the data points from each of the three distributions 
collapse onto a single curve, confirming equationj^ This result has also been confirmed 
using experimental size spectra from a cirrus cloud — see Westbrook et al (2004a,b) for 
details. 

In summary, the crucial results (backed up by experimental data) are that the rela- 
tionship between mass and linear dimension takes a power law scaling, with a universal 
exponent determined by the physics of the aggregation process, and a non-universal 
prefactor which is sensitive to both the pristine crystal type and size. In addition, the 
aggregate size distribution is described by a universal underlying distribution function, 
which is simply rescaled as a function of the average particle mass as the aggregation 
proceeds. 



4. Radar cross sections oe individual aggregates 

{a) The Rayleigh-Gans theory 

In this section the back scattered intensity from a single snowflake is calculated 
using the Rayleigh-Gans theory. The essence of this approximation is that the particle 
is split up into a number of small volume elements dv at position r relative to the 
overall centre of mass. Each element is treated as a Rayleigh scatterer, and we ignore 
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any interactions between the elements (each one sees only the applied incident wave). 
Summing up the contributions from each of them with an appropriate phase factor, the 
radar cross section of the complete particle is obtained (see Bohren and Huffman 1983): 
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Here pi^e and e are the density and dielectric constant of (solid) ice respectively, e- 
the direction of propagation of the incident wave, and k = 2t: /\. 

The form factor / represents the deviation from the Rayleigh regime as the overall 
size of the particle and the wavelength of the incident light become comparable. In 
the Rayleigh limit where the particle is much smaller than the wavelength, / = 1, and 
the radar cross section ^ takes the well known oc w?X~'^ dependence. As the particle 
size increases relative to the wavelength, the form factor falls off with a functional 
dependance which, in general, depends on the geometry of the particle whose volume is 
being integrated over in equation|6l 

The key advantage of the Rayleigh-Gans approach is it allows the scattering cross 
section for a particle of any shape to be calculated, provided the particle is not too 
large or too strong a dielectric. Unlike the discrete dipole approximation (Purcell 
and Pennypacker 1973; Draine and Flatau 1994), the coupling between the elements 
is neglected, and this places constraints on its applicability. In particular, Bohren 
and Huffman (1983) quote the conditions \ y/e — 1| < 1, \y/e — l\kr < 1. However, for 
fractal aggregates (such as our snowflakes). Berry and Percival (1986) have shown that 
these conditions may be significantly relaxed. Their calculations show that because 
of the open structure, multiple scattering between the monomer particles composing 
fractal aggregates with dj ^2 is negligible, irrespective of the overall size of the 
aggregate, and for > 2 only becomes significant when the number monomers per 

aggregate becomes of order (fcro)^'^^/*^'^^"^''. Thus provided that multiple scattering 
within the monomers themselves is negligible (ie. \^/e — l\krQ < 1), which we expect 
to be reasonable for typical radar frequencies/monomer sizes, the Rayleigh-Gans theory 
may be safely employed. 

The pristine crystal geometry may also play a role in determining the range of appli- 
cability of the Rayleigh-Gans theory. Since all interaction between the scatterers (which 
are assumed to act as equivalent volume Rayleigh spheres) is ignored, the anisotropy of 
the monomer particles is implicitly neglected. However, Liu and Illingworth (1997) have 
studied the scattering from a quite elongated hexagonal ice crystal (aspect ratio=3), and 
found that for kr^ ^ 1 an equivalent volume sphere gives the same radar cross section 
to within 1%. We therefore expect that ignoring the anisotropy of the ice crystals is an 
acceptable approximation, provided that they are small enough compared to the incident 
radar wavelength. 



(b) The form factor 

It is apparent from equations |5] and |6l that outside the Rayleigh limit the radar 
cross section is sensitive to both particle size and shape through the form factor /. For 
every aggregate produced in our simulations, equation Q was evaluated for a range of 
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Figure 4. The form factor / plotted as a function of 2kr for our synthetic aggregates. Left hand panel 
shows the form factor evaluated for all of the aggregates produced over the course of a simulation, binned and 
averaged as a function of 2kr (solid line with dots). The dashed line shows the small size Guinier approximation 
/ ~ 1 — i{2fer)^; the dotted line shows the form factor for a sphere with the same radius of gyration r (note 
r ~ 0.3D). The right hand panel highlights the insensitivity of / to pristine particle shape, as the aspect ratio of 
the pristine columns is varied from unity (solid line) to 0.25 (dashed line) with no appreciable change in the curve. 



wavelengths. The resulting values of / were binned as a function of {2kr) and averaged 
to yield the curve in figure 0^. 

Given our results on the universality of aggregate snowflake geometry, we anticipate 
that the form factor will also be independent of the details of the monomer crystals: in 
figure we illustrate this point by varying the aspect ratio of the pristine crystals and 
measuring an almost identical form factor. 

We now compare the form factor of our aggregate snowflakes with analytical results 
from the literature. It is well known (particularly in the polymer physics community) that 
the asymptotic departure from the Rayleigh limit at small sizes may be approximated 
as: 

/~l_i(2fcr)2. (7) 

This is the Guinier equation (Guinier 1939, Guinier et al 1955), and it applies to particles 
of any shape, provided 2kr is small enough. Equation is plotted alongside the 
aggregate form factor in figure |^ and proves to be an excellent fit to our simulation 
data up to 2kr ~ 1, beyond which it underestimates the aggregate curve. 

In meteorological studies, ice particles are often approximated as spheres, in order 
to utilise Mie's exact solution for the scattering from a dielectric sphere (Mie 1908). The 
form factor for a sphere is also well known (Bohren and Huffman 1983) and is given 
by: 

2 

(8) 

where u = 2\/b/2>kr. Matrosov (1992) has employed this formula as an alternative to 
the more complex Mie expansion, and showed the results closely mimic the full Mie 
solution. 

Figure |4] shows that the correspondance between the spherical model and our 
aggregates is quite good at small enough sizes, but diverges once 2kr > 0.5. We also 
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note that the sphere model predicts a form factor of zero (ie. no back scatter at all) at 
2kr ~ 3.5, and similar points with / = at higher values of 2kr. This undulating form 
(also found in the exact Mie solution) is a peculiarity of the spherical case, and is not 
reproduced for our aggregates which do not possess the smooth, symmetrical shape of 
the idealised dielectric sphere. The errors associated with the spherical approximation 
are discussed in more detail in section 8. 

Values of 2kr for typical Cirrus snowflake sizes (0.5 and 1.0mm maximum dimen- 
sion) at representative radar frequencies (35 and 94 Ghz) are marked on figure |4j and 
in all cases 2kr < 1. We investigate the consequences of this in the next section, where 
we show that if the Guinier approximation holds, derivation of average particle size is 
exceptionally straightforward and may be achieved through a simple analytical formula. 

Since radar scattering is rather sensitive to the largest particles in the distribution 
(even though there may be relatively few of them), we also consider the case of 
snowflakes beyond the Guinier regime and fit a curve to the complete form factor. We 
then demonstrate how such a curve may be used to interpret radar data which is outside 
the Guinier regime. 



5. INTERPRETATION OF REFLECTIVITY DATA 

In this section we employ the dynamical scaling property of the distribution, along 
with the results from the previous section on the scattering from individual snowflakes, 
in order to study the physical significance of reflectivity measurements and to make use 
of this understanding to infer the cloud's microphysical properties (average snowflake 
size, ice water content, etc). 

Having calculated the form factor /, equation |5] may be evaluated for each 
snowflake in the scattering volume and the results added, to obtain the reflectivity 

^=(^)E"^'^(2A:r) (9) 

where the constant c = 367r^|^p/p|^g. Multiple scattering between the snowflakes 
themselves is assumed to be negligible, given the dilute nature of most ice clouds. 
Note that the reflectivity is commonly rescaled to give the 'radar reflectivity', defined 
as Z = A^7r~^|^:=^|^^?7. In what follows we will use rj for clarity, since its definition 
follows naturally from the radar cross sections calculated in the previous section. 

(a) The Guinier regime (small aggregates) 

From equation the information contained in measurements of the radar reflec- 
tivity is immediately apparent. If the snowflakes are much smaller than the wavelength 
then / ~ 1 and we recover the well known Rayleigh limit t] oc A^^ Y1 ^^id the re- 
flectivity is simply a measure of the second moment of the mass distribution. As the size 
of the snowflakes and the wavelength become comparable, the particle radius also con- 
tributes to the reflectivity through /. Provided the combination of snowflakes size/radar 
wavelength falls within the Guinier regime 2kr < 1, the form factor is approximated by 
equationQand as a result the reflectivity is simply: 

The first term in the bracket is the Rayleigh result, and for long enough wavelengths 
(compared to r) this behaviour dominates the scattering. As wavelength and particle 
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size become more comparable the second term becomes more significant, and the 
scattering becomes directly dependent on the particle radius. The implication is that 
if more than one radar wavelength is used, an estimate of average particle radius may 
be obtained. Given equation ^| a natural definition for that average radius is Vav = 
[Y^ rr^r^ I ^ m^]2, and for two reflectivity measurements r/i and r]2 (at wavelengths 
Ai and A2) we find: 



1-/3 



(11) 



167r2 1/A2-/3/A; 

where /? = (r?iAf ) lij]2^^, or in terms of radar reflectivities /5 = Z1/Z2. 

The practical consequence of equationfTDis that given two simultaneous measure- 
ments at different radar frequencies the average snowflake radius may be calculated 
analytically through equation^^ provided that a) at least one of the radars is operating 
at a wavelength sufficiently short that some of the snowflakes fall outside the Rayleigh 
regime, and b) both radars have wavelengths long enough that the snowflakes fall pre- 
dominantly inside the Guinier regime (2fcr < 1). In previous dual wavelength studies 
(eg. Matrosov 1998; Hogan et al 2000), the relationship between average particle size 
and /? has been calculated numerically for spheres and ellipsoids; here we have a simple 
analytic result, which we expect to be applicable to typical snowflake sizes and radar 
wavelengths. 

Once Tav has been calculated then either of the reflectivity measurements may be 
used to infer ^ rr? . Rearranging eguationfTUIvields: 



m2 = !A ^2) 

c{\-\{2kravf\ 



In general we would like to use Vav and ^ to estimate the ice water content 
IWC = ^ m. From our definition of the average mass we see that this is simply given 
by: 

IWC = ^^ . (13) 

Thus the IWC may be inferred, provided a measure of average snowflake mass is 
available. Unfortunately it was demonstrated above that dual wavelength measurements 
only give a measurement of average radius, not average mass. However equations|3]and 
Elprovide a means to covert between the two: 

■mav = a— ^'L (14) 

(taking df = 2). The ratio P2/P3 has a universal (dimensionless) value, which we 
measure in our simulations to be around 0.43 ± 0.01; analysis of experimental size 
distribution data (see appendix B) agrees closely, with a value of 0.44 it 0.01. The 
value of the prefactor a by contrast is not universal, depending upon both the size and 
geometry of the pristine ice crystals making up the snowflakes. For simple, compact 
crystal habits such as columns or plates there is a linear sensitivity to the crystal size 
a oc ro; for dendritic shapes, which have a fractal scaling of their own, we expect the 
prefactor to be less sensitive, since it is the mismatch between the mass-radius scaling 
of the crystals and that of the overall aggregate which results in the dependence on 
crystal size. 
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(b) Larger aggregates 

For large enough snowflakes relative to radar wavelength, the Guinier formula 
breaks down, and in this section we identify the point where this occurs, and attempt 
to fit a curve to describe the scattering at large sizes. In order to achieve this, we rewrite 
equation|9lin the form: 

^ = ^ (E"^') ^(2A;r,,) (15) 

defining 9"= "i^// Yl (i^- m^-weighted average of the form factor /). We 
have already identified the form factor for individual aggregates as being a universal 
function of 2kr. Because of the dynamical scaling of the distribution and the universal 
power law scaling between mass and radius, we expect 9" be a universal function of 
2krav Expressing r] in this way makes it immediately obvious that only two quantities 
may be directly derived from the reflectivity: ^ and rav AH other quantities (eg. 
ruav, IWC) must be inferred from these two, and any other available measurements or 
assumptions. 

Figure |5] shows 3' calculated from our simulations. In the Guinier regime, the 
average form factor is simply 1 — i(2fcra„)^, and this approximation fits the 
simulation data well up to 2krav — 1, beyond which it underestimates the scattering. 
To try and provide a description for the scattering beyond 2krav ^ 1 we have fitted a 
curve to our simulation data. At small sizes we know that the curve must approach the 
Guinier formula; at large sizes, it is well known in the physics literature (eg. Viscek 
1989) that Joe (2A:ra„)-'^^ We therefore fit a curve of the form: 

J = l + ci(2fcr,.)2 

l + (ci + i)(2A:r,„)2 + c2(2fer,,)4 

which has the correct asymptotics in both limits. As shown in figure |5l this provides a 
good fit to the simulation data with ci = 12.7 and C2 = 3.6. 

Given a dual wavelength ratio measurement (3, we may use this fitted curve to 
estimate the average radius. Noting that j3 = 3^{2kirav)/3^{2k2rav) (where ki = 27r/Ai, 
k2 = 27: /X2) we obtain: 

a'i4rlf + 6'(4rL)2 + c'(4rL) +d' = (17) 

where: 

a' = ciC2klkl{(5kl - kl) (18) 

b' = C2iPkt-ki) (19) 

c' = ci{pkj-kl)+ (20) 

+ (ci + i) (/3A;2 - A:2) (21) 

d' = 13-1. (22) 

EquationlTTlis a cubic equation in 4r^^, and may be solved using standard methods (eg. 
Press et al 1992) to obtain Vav The second moment may then be obtained straightfor- 
wardly, using either of the reflectivity measurements: ^ m? = rjX^ /c3^{2krav)- Inter- 
pretation of these results in terms of IWC then follows as described in section (a). 



6. Inferring other microphysical properties 



Moments of the size distribution other than the ice water content may be also 
be inferred from reflectivity measurements. In section 5 it was demonstrated that 
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Figure 5. The m-^ -weighted average form factor for the whole distribution 3^=Y1 ™-^// Yl plotted as 
a function of average size 2krav for our synthetic aggregates (points w/solid line). Dotted line is the Guinier 
approximation 1 — ■|(2A,Tai,)^; dashed line is our fitted curve (see text for details). 



dual wavelength measurements of ice clouds allow the second moment to be 

measured, along with an average radius Vav, which may be converted to an average mass 
niav given knowledge of the prefactor a in the fractal scaling relation. We therefore wish 
to relate the moments of the size distribution that are of interest to the second moment 
(which we know). From the dynamical scaling of the distribution @ it follows that: 



m 



m„ 



-2P' 



P2 ^ 



m 



(23) 



forn ^ 1. Moments with n < 1 scale differently due to the power law at the small end of 
the distribution — see Westbrook et al 2004b for details). The ratios Pn/P2 are universal, 
and may be measured from simulations or derived from experimental size spectra as per 
Appendix B. 

A parameter of significant interest to meteorologists which some authors have 
sought to derive from radar measurements is the effective radius rg. This characterises 
the average ratio of particle volume to projected area (Foot 1988): 



3IWC/P,, 



(24) 



where A is the particle projected area. For aggregates however, estimating from the 
dual wavelength derived average size is not possible. Since the mass is proportional to 
r^, and likewise for the projected area (the aggregates considered here are fractals with 
df = 2), the effective radius does not scale with the radius of gyration or maximum 
dimension; its value is governed by the details of the pristine particles that compose 
the aggregate snowflakes. Although there is likely to be some correlation between r^^ 
and Te from cloud to cloud, this simply reflects the fact that larger pristine particles 
tend to yield larger aggregates. Essentially rav is a measure of the overall scale of the 
aggregate, while r^. is a measure of the size of the pristine crystals. Trying to derive one 
radius from the other is therefore not possible, and dual wavelength radar measurements 
do not provide the means to estimate r^. 
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Figure 6. Time series showing the dual wavelength ratio /3 = •Z35Ghz/^94Ghz between 15:00 and 19:30 on the 
22"<^ of June 1996 measured over Chilbolton in the UK (top panel). Shown underneath is the inferred average 
radius rav derived using both the analytic Guinier formula (dashed line) and the complete fitted curve of section 

5b (solid line). Units are in millimetres. 



7. Application to cirrus cloud data set 

In this section we apply the two methods for estimating rav to a dual wavelength 
dataset from a real ice cloud. On the 22"^^^ of June 1996, vertically pointing radar 
measurements of a cirrus cloud over Chilbolton in the UK were made using 35 and 
94Ghz radars. This data was presented in Hogan et al (2000), and the full experimental 
details are given in that paper. Here we analyse the radar reflectivity time series from 
one representative altitude (~ 6 km) where the data shows approximately the full range 
of reflectivity values measured during the experiment. From these measurements we 
have calculated the dual wavelength ratio j3 = Z^^^qi^^ / Zgichz and this is plotted as a 
function of time in figure |6l From this we have used first the Guinier approximation 
(small flakes, as described in section 5a) and then our fitted curve (section 5b) to infer 
the average radius Vav These results are also plotted in figure |6l and it is immediately 
apparent that there is very little to distinguish the Guinier-derived radius from the one 
obtained using the full fitted curve. The implication is that the snowflakes are small 
enough that they lie within the Guinier regime, and since the average radius Vav only 
reaches around 0.3mm at most (~ equivalent to 1mm maximum dimension), this means 
'^krav ^1.2 for 94Ghz, and 2krav ^ 0.4 for 35Ghz. Refening to section 5 these values 
both fall roughly within the Guinier regime. 
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Figure 7. Dual wavelength ratio /3 as a function o f av erage radius Vav in millimetres for 35- and 94-GHz radars. 
Black line shows the prediction of our fitted curve 116t ; dashed line shows the Guinier approximation. Left hand 
panel shows good agreement between the Guinier approximation and our fitted curve at small sizes (r < 0.25); 
the right hand panel shows the behaviour at larger sizes, where there is a marked divergence. 



For the studied cirrus cloud at least, the Guinier formula is a good approximation. 
We note however that even here the largest particles were on the edge of the regime 
{2krav — 1); for thicker clouds containing larger aggregates this formula may not be 
sufficient, and the complete fitted curve ( fT6t should be used. Figure shows a plot of 
dual wavelength ratio as a function of average radius for 35- and 94-Ghz radars, as 
predicted using both (fT6t and the Guinier approximation. It is apparent that the Guinier 
formula holds up well to around Vav — 0.3mm, beyond which point it rapidly breaks 
down. It is useful to see however, that at least for parts of the cloud dominated by smaller 
aggregates, the Guinier approximation does hold and that a simple analytic expression 
provides an accurate estimate for In both cases there is a need for the methods 
to be validated through simultaneous aircraft measurement of particle mass and radius 
distributions. 

8. Comparison WITH SPHERICAL MODELS 

Many authors have used a spherical approximation when estimating the back 
scatter from snowflakes, modelling them as homogeneous mixtures of air and ice. The 
motivation for this kind of approximation is that the scattering from a sphere with a 
given diameter and dielectric constant may be calculated exactly using Mie theory. Here 
we consider how closely the spherical model matches the results from our synthetic 
aggregates. 

To mimic the results for the radar cross section by our simulated aggregates using 
a sphere, we attempt to match up both the form factor / for the aggregates, and the 
dependance for 2A;r ^ 1. In figure |4] the form factor for a sphere was plotted alongside 
that for our simulated aggregates. There is good agreement up to 2kr ~ 0.5, after which 
the curves diverge somewhat. Thus for small aggregates at least, an equivalent sphere is 
an accurate representation provided that the mass and radius of gyration are conserved. 
We focus therefore on this regime, and consider the errors at larger sizes later. In order to 
match up the scattering from our sphere with that of the aggregate it is intended to model, 
we must ensure the radius of gyration of the sphere is the same as that of the aggregate. 
For our simulated aggregates, the radius of gyration and maximum span are found to 
scale linearly with r ~ 0.3D, whereas a sphere of diameter Dgp has r = 0.39Dsp (see 
Appendix A). To match these up then requires that the diameter of the equivalent sphere 
be Dsp = 0.77D. This ensures that the sphere has the same radius of gyration as the 
aggregate it is intended to represent, and hence the form factor shown in figure |3 A 
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Figure 8. Dual wavelength ratio /3 as a fu nctio n of rav for 35- and 94-GHz radars as per figure Black 
line shows the prediction of our fitted curve <16> : dashed is the Guinier approximation. The grey lines show 
measurements of 13 from a computer simulation in which a sphere was used to model the aggregates, where: (i) 
the radius of gyration of the aggregate and the sphere used to model it were matched (light grey), and (ii) the 
diameter of the sphere was chosen to be the maximum dimension of the aggregate (dark grey). 



different choice of Dsp would result in the curve being squashed or stretched along the 
2kr axis, and under/over-estimating the back scatter for 2kr < 0.5. 

To estimate the error involved in this approach we have run computer simulations 
where we calculate the dual wavelength ratio at 35- and 94-Ghz, modelling the scattering 
from our aggregates as that of spheres. We have considered two cases: (i) where we 
choose each sphere's diameter so as to match the radius of gyration of the sphere with 
that of the aggregate as discussed above, and (ii) where we choose the diameter of the 
sphere to match the maximum dimension of the aggregate D (a model suggested by 
some authors in the literature): the results of these simulations are shown in figure |H] 
When the radii of gyration are matched, the agreement at small sizes is good, with only 
around a 5% error in (3 relative to our fitted curve at rav = 0.25mm. At larger sizes the 
divergence is more apparent, and the sphere model overestimates (3 by around 25% at 
Vav = 1mm. Model (ii) performs rather less well: matching the diameter of the sphere to 
the maximum dimension, the dual wavelength ratio is overestimated by around 20% at 
fav = 0.25mm, and by more than 100% at Vav = 1mm. We conclude that both methods 
overestimate f3 for a given Vav, but matching the radius of gyration of the sphere with 
the aggregate results in much lower errors than matching their maximum dimension. 

The second dependence that must be matched is the snowfiake mass. Since the 
snowfiakes are not solid ice spheres, information about how much of the mixture is air 
and how much is ice must be introduced by using an 'average' dielectric constant e^^ 
weighted by the volume fraction of ice in the sphere fy = m / ^{D spY Pice- The most 
common prescription for €av is the Maxwell-Garnett (1904) formula: 

_ 1 - A + <fv 

" 1 - /. + a. ^ ^ 

where the parameter C, reflects the distribution of sizes and shapes of the ice 'inclusions' 
in the mixture (see Bohren and Huffman 1983 for more details). For spherical inclusions 
C = 0.58. To check that introducing the aggregate density in this way yields the same 
m? dependance as per equation |5j we study the Rayleigh limit where the back scatter 
cross section of our equivalent sphere is given by: 
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Figure 9. Ratio comparing the back scatter cross section calculated using the Maxwell-Gamett sphere model, 
and the Rayleigh-Gans theory, for 2kr <^ 1 (the Rayleigh limit). Solid line represents spherical inclusions 
(k = 0.58), dashed is Meneghini and Liao (k = 0.63), and dotted is rod shaped inclusions (k = 0.65). 



For solid ice spheres eav = e, and since m = ^{DgpYpice, one recovers equation |5] with 
/ = 1. For small volume fractions (the case most relevant to fractal aggregates, since 
the overall density is proportional to r^^~^) we calculate the ratio of the back scatter 
obtained from equation|26lto that obtained from equation|5las a function of the volume 
fraction fy. in the Rayleigh limit this turns out to be simply | ^°''^2 l^/{|f^P/^}- This 
ratio is unity at = 1; even as /„ — > the error is only around 1% (see figure |9ll. 
Other values of the parameter C, are also in common usage - the assumption of needle 
shape inclusions gives C, = 0.65 and Meneghini and Liao (1996) have proposed a value 
oi C,= 0.63. These alternate values of C, essentially represent the polarisability of the 
inclusions (which we interpret in the case of our aggregates as being the monomer 
particles), and they predict a back scatter increased by as much as 25% as fy — > 0, as 
shown in figure|9l In section 4 it was commented that in the Rayleigh limit the scattering 
from the individual monomers should be well represented by a sphere of the same 
volume; if this is the case then the Rayleigh-Gans approximation ought to be accurate, 
and these two values of Q overestimate the scattering quite significantly, particularly 
at low volume fractions. This is perhaps a reflection of the fact that 'effective medium' 
theories such as the Maxwell-Garnett rule are intended for strong dielectrics at relatively 
large volume fractions, neither of which is really the case for ice aggregates. 

In summary, the spherical approximation matches the results for our aggregate 
snowflakes provided that i) the diameter is chosen so as to match the radius of gyration 
of the sphere to the snowfiake it is intended to model, and ii) the parameter C, is chosen so 
as to give the correct m?' dependence in the Rayleigh limit. If the size of the aggregates 
is not too close to the wavelength of the radar, then the error in /3 is only 5% or so; for 
larger aggregates this may increase to 25% or more. It is important to note that a 5% 
error in (3 may translate into a rather more significant error in the derived parameters 
such as r„„ and IWC. 




RADAR SCATTERING BY AGGREGATE SNOWFLAKES 



15 



9. Conclusions 

The radar cross sections for our simulated snowflakes have been calculated using 
the Rayleigh-Gans theory by treating the aggregates as being an assembly of inde- 
pendant Rayleigh scatterers, and summing up the contributions from each element in 
that assembly. In this theory, the scattering is proportional to the square of the particle 
mass (as in the Rayleigh limit) and the form factor /, a dimensionless function which 
describes the depedence on particle size and geometry. This function is universal (ie. 
independent of the size and shape of the monomer particles), and / has been calculated 
for our simulated aggregates. 

This is (to the authors' knowledge) the first time that the scattering from realistic 
aggregate snowflakes has been calclulated. It also appears to be the first time the 
Rayleigh-Gans theory has been employed to study radar scattering by ice particles, 
with the exception of Matrosov (1992) who used it only as an approximation to the 
exact Mie solution for a sphere. Ideally we would like to verify the results obtained in 
this paper using the discrete dipole approximation which accounts for the interaction 
between the scattering elements. However, this places heavy demands on computer 
time and memory, since a large number of dipoles are needed to accurately model the 
detailed structure of the aggregates, and calulations for many possible realisations and 
orientations must be done in order to obtain good statistics. So far this has not been 
achieved. Also, in order to test the accuracy of the theoretical results in sections four 
and five, simultaneous radar and aircraft observations of ice clouds are needed. 

The form factor for our simulated snowflakes was compared to analytical results for 
/ from the scattering literature, the most successful of which in describing the data was 
the Guinier result / = 1 — ■|(2A;r)^, which fits our data well up to 2kr ~ 1. The result 
for the sphere is a good approximation below 2kr ~ 0.5, beyond the curves diverge 
somewhat, the sphere first over- then under-estimating the back scatter. 

Having calculated the scattering cross section of individual aggregates, these were 
then integrated over the size distribution, and the results interpreted in terms of moments 
of that distribution. The reflectivity is proportional to the second moment of the mass 
distribution ^ and the average form factor 9^ = J2 f / J2 . We have calculated 
as a function of the average radius Vav, and fitted a curve to our simulation data. The 
application of our results to the inference of average snowflake size from radar data has 
been investigated, and for snowflakes sufficiently small compared to the wavelength, the 
Guinier approximation provides a very simple analytic expression for Vav in terms of the 
dual wavelength ratio /3. At larger values of 2krav, the Guinier curve underestimates the 
back scatter: we therefore fitted a curve to our simulation data and used that curve to 
provide a method to interpet the dual wavelength data. 

Ideally, we would like to interpret the radar reflectances in terms of the ice water 
content (total mass per unit volume) in the cloud. Since reflectivity is proportional 
to ^ m^, some measure of average particle mass is required. We have shown that it 
is possible to derive the average particle radius Vav from dual wavelength radar data 
and this may be converted to an average mass, provided that some measurement or 
prescription for the prefactor a in the fractal scaling relation is available: a is a function 
of the pristine crystal geometry and size. 

Some authors have attempted to derive the effective radius (characterising the ratio 
of particle volume to projected area) from radar measurements. However we have shown 
that Te depends only on the monomer particles (to which the radar is insensitive) and 
not the overall size of the aggregate as characterised by rav Direct inferral of Vg then is 
not possible from dual wavelength radar data. 
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Finally, the errors involved in modelling our simulated snowflakes as air/ice spheres 
using the Maxwell-Garnett mixture theory were analysed. It was found that provided 
some care was taken in constructing the equivalent sphere, and the radius is not too large, 
the error is around 5% {Vav = 0.25mm with 35/94Ghz radars), increasing to around 25% 
at larger sizes (r^^ = 1mm). It is worth noting that such an error in the dual wavelength 
ratio has the potential to translate into a much larger error in the derived parameters, 
especially if f5 is close to unity. 
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Appendix A 



The radius of gyration 

In this paper we have used the radius of gyration r as the characteristic aggregate 
length scale. It is defined by splitting the particle into small volume elements, with 
position r and mass dm. Then: 



f |rpdm 
/ dm 



(A.l) 



integrating over the particle volume. This radius is a natural length scale to use for the 
scattering calculations because it is closely related to the integral for / in equation 6 
at small kr, leading to the Guinier expansion / ~ 1 — i(2A;r)^. For our aggregates it 
is linearly related to most other characteristic length scales, in particular the maximum 
dimension D of the snowflakes (r ~ 0.3D). 



Appendix B 

Inferring the universal moment ratios Pn/P2ffom experimental snowflake span 

distributions 

Here we show how the universal ratios P2/P3 (used to link mav to Vav) and Pn/P2 
(used to convert ^ and mav into a different moment of the distribution "^"') 
may be calculated from experimental snowflake span distributions. Using equation |3l 
the moments of the radius distribution M„ are given by: 

/•oo 

Mn= {m/a)''/'^'n{m,t)dm (B.l) 
Jo 

which, using the dynamical scaling property @, and taking = 2 is: 

M„ = (m,,/a)"/2-€p„/2 (B.2) 

From the definition of mav, one has P2/P1 = 1- We therefore obtain the ratios M4/M2 = 
mav/cL, and M6/M4 = {mav/a){p'i/p2)- The ratio we seek to calculate, P2lpz, is then 
simply: 

P2 _ Ml 
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We note that it makes not difference whether we use moments of the radius distribution 
(r) or of the span distribution (D) since the end result is dimensionless, and r (x D. 
We have calculated P2/P3 from the experimental particle span distributions presented in 

Westbrook et al (2004a,b), and the ratio is found to be 0.44 ± 0.02. 

The ratio Pn/P2 may be obtained in a similar fashion. From equation IB . 21 above : 

M2n (mav\-^-'^ Pn 

(B.4j 



M4 \ a J p2 
and since M4/M2 = niav/a, the ratio we want is simply: 



1 • (B.5) 

P2 M^-^ 
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